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We estimate the photon and dilepton emission rates from 
hot hadronic matter with in-medium spectral shift and broad- 
ening of vector mesons. It is observed that both the WA98 
photon data and CERES /NA45 dilepton data can be well re- 
produced with similar initial conditions. The freeze-out con- 
dition has been constrained by the transverse mass spectra of 
pions and protons measured by the NA49 collaboration. We 
argue that simultaneous measurement of the pr spectra of 
single photons as well as invariant mass distribution of dilep- 
tons is crucial to understand the in-medium spectral function 
of the vector mesons. 

PACS: 25.75.+r;12.40.Yx;21.65.+f;13.85.Qk 



1. INTRODUCTION 

The study of the behaviour of hadrons at finite tem- 
perature and/or high baryon density has attracted sub- 
stantial theoretical and experimental attention in recent 
times. The theoretical motivation stems from the ex- 
pectation that broken chiral symmetry of QCD might be 
restored at high temperature and/or density, at least par- 
tially and this has nontrivial effects on the spectra of low 
lying hadronic states. Again, these studies are of consid- 
erable importance in connection with the experimental 
detection and study of quark gluon plasma (QGP) in 
the ultra-relativistic collisions of heavy ions. This is be- 
cause disentangling the signals from quark matter would 
require a precise estimate of emissions from hadronic 
sources which constitutes an overwhelming background 
in such experiments. Now, the properties of short-lived 
resonances like the p and u> mesons which are likely to de- 
cay within the fireball can only be studied through deep 
probes such as photons and lepton pairs which are essen- 
tially electromagnetically interacting and thus probe the 
entire space-time volume of the collision [1]. The medium 
modifications of low mass vector mesons has been the 
subject of numerous theoretical investigations culminat- 
ing broadly into the following two types of predictions. 
QCD based models predict a downward shift in the mass, 
which we will mention as spectral red-shift and a rescat- 
tcring scenario giving rise to an enhanced width which 
we have denoted by spectral broadening, the mass re- 
maining largely unchanged. Our aim in this work is to 
calculate the electromagnetic spectra considering various 
plausible scenarios and compare with the WA98 [2] pho- 
ton and CERES/NA45 [3] dilepton spectra in order to 
comment on the relative importance of the two scenarios 



mentioned above. 

Our analysis will proceed as follows. In the next sec- 
tion we will define the effective mass and width that have 
been used in the calculation and define the scenarios of 
medium effects considered. At the onset we emphasize 
that our main objective is to see the effect of these con- 
trasting scenarios on the electromagnetic spectra and not 
on the details of how the forms of the effective mass and 
the width have been arrived at. Consequently we will 
use phenomenological arguments to justify our choice of 
the parametrisations for these quantities. In section 3 
we will discuss the initial conditions in considerable de- 
tail. In the relevant transverse momentum region of the 
photon spectra the main sources of photons are the hard 
pQCD processes as well as the thermal emissions. Con- 
sequently, the number of participants as well as the ini- 
tial temperature are the relevant initial conditions which 
have been estimated. For the low mass dilepton spectra, 
the contribution from Dalitz decays as well as thermal 
emissions are important. Since photons and dileptons 
are emitted at all stages of the collision, the space-time 
evolution plays a very significant role in estimating the 
total yield. We have used two different approaches to 
study this aspect; relativistic hydrodynamics and trans- 
port. As we shall see, medium effects do play a very 
important role in the equation of state and consequently 
the cooling profile. These aspects have been dealt with 
in section 4. The transverse mass (mx) spectra of pi- 
ons and protons obtained at SPS have been evaluated 
and compared with data from NA49 [4] in section 5 in 
order to estimate the freeze-out condition which goes as 
an input into the hydrodynamic calculations. With all 
these inputs we evaluate the photon and dilepton spec- 
tra in sections 6 and 7 respectively, section 8 contains a 
summary and discussions. 



2. EFFECTIVE MASS AND WIDTH OF 
HADRONS 

It has been emphasized that the properties of hadrons 
will be modified due to its interaction with the particles 
in the thermal bath and such modifications will be re- 
flected in the dilepton and photon spectra emitted from 
the system (see Refs. [5-9] for review). Broadly two types 
of medium modifications are expected: shift in the pole 
position and/or broadening of the spectral function. As 
discussed earlier, our main objective is to comment on 
the nature of medium effects through the study of WA98 
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photon and CERES dilepton data. Accordingly we con- 
sider the following scenarios: (I) the system is formed 
in the hadronic phase with the hadronic masses (except 
pseudoscalars) approaching zero near the critical temper- 
ature according to the universal scaling law [5,10]. The 
in-medium masses of the hadrons and the decay width of 
p are taken as 

m* H /m H = (l-T 2 /T 2 f 

m*Jl - 4m 2 /mff' 2 (1 + 2f BE ) (1) 
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where /be is the Bose-Einstein distribution. For the 
width of lj see Ref. [9] . (II) the width of the vector meson 
(p) increasing with temperature as, 



p 



V (1 - T 2 /T 2 ), 



(2) 



and the masses remain constant at their vacuum values. 
The value of A in eq. 1 is 1/6 and 1 /2 for Brown- Rho and 
Nambu scaling [5] respectively. The present experimen- 
tal statistics can not differentiate among different values 
of A. Therefore, we leave it as a parameter here. The 
scenario II derives motivation from the results of chiral 
models according to which the mass of the p does not 
change to order T 2 but the leading order T-dependcncc 
of the pion decay constant results in the p decay width 
having the form described by Eq. 2. We have also con- 
sidered a third case, (III) in which both the masses and 
widths are maintained at the vacuum values. In eqs. 1 
and 2, T c is the critical temperature of chiral symmetry 
restoration where the quark condensate (qq) and conse- 
quently the mass of the hadrons tend to zero. 

Let us study the in-medium effects on the thermal 
phase space factor through the p meson as it plays the 
most important role for the electromagnetic probes. In 
Fig. 1 we display the change in the density of thermal p 
as a function of temperature for scenarios (I), (II) and 
(III) by solid, dotted and long-dashed lines respectively. 
It is clear that the density of p is very sensitive to the 
spectral shift of p mass due to Boltzmann enhancement 
but rather insensitive to the change of width. The reason 
for small change in the density of p due to the in- medium 
broadening can be understood from the following argu- 
ments. The density of an unstable vector meson (say, p) 
in a thermal bath can be written as [11], 
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<Pkd?xds (2tt) 3 e VW+7s/T _ i 



P{s) 



(3) 



where g is the statistical degeneracy of the particle and 
P(s) is the spectral function, 



P(s) 



1 



imn 



7r(s-TO 2 -RcII) 2 + (Imn) 2 



(4) 



Imn (Ren) is the imaginary (real) part of the p self 
energy. Eqs. 3 and 4 indicate that the density of par- 
ticles (vector mesons) in a thermal bath is given by the 



Bose-Einstein distribution weighted by the Breit-Wigner 
function, which gets maximum weight from the value of 
s = m 2 + Re n, the contribution from either side of the 
maximum being averaged out. Therefore, the results be- 
come sensitive to the value of s — m 2 * = m 2 + RelT and 
not to he width of the spectral distribution. Note that 
P(s) tends to 5(s — to 2 — RcII) as Imll — > 0, corre- 
sponding to a stable particle. (Here, II is proportional to 
the trace of the self energy tensor U pL/ of the p) . 
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FIG. 1. Solid (long dashed) line represent density of ther- 
mal p as a function of temperature with in-medium mass and 
vacuum width (vacuum mass and width). Dotted line repre- 
sents the density with vacuum mass and in-medium width. 



3. INITIAL CONDITIONS 

Here we will discuss the parameters which enter as in- 
puts to the evaluation of electromagnetic and hadronic 
spectra from relativistic heavy ion collisions correspond- 
ing to different scenarios considered. We begin with the 
quantity n(b) which is the average number of nuclcon- 
nucleon collisions in the collision of two nuclei of mass 
numbers A and B at an impact parameter b. This is 
given by [12], 



n(b) = ABT AB (b)a v 



(5) 



where Tab (b) is the nuclear overlap integral evaluated 
from the following expression, 



T AB (b)= J d 2 sT A (s)T B (b-s), 



(6) 



and the thickness function, T A is defined as T A (b) = 
Jdzp A (b,z). The nuclear density, p A (b, z) is 
parametrized by Wood-Saxon type profile function [13], 



Pa(t) = Pa 



1+ujt 2 /R 2 a 

1 _|_ e (r-R A )/z 



(7) 



where u> — 0, z = 0.549 for the Pb nucleus and R A = 
\.2A X I Z and the central density, po is determined by the 
normalization condition J d 3 rp A (r) = 1. In Fig. 2 n(b) 
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is shown as a function of the impact parameter, b. The 
most central event in WA98 corresponds to b ~ 3.2 fm/c. 
At this value of the impact parameter n(b) <~ 660. 
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FIG. 2. Hadron multiplicity, effective number of nu- 
cleon-nucleon collisions and initial temperature as a function 
of impact parameter calculated by using Glauber model for 
nucleus - nucleus collisions (see text). 

The hadron multiplicity in A + B collisions at an im- 
pact parameter b can then be evaluated in terms of the 
multiplicity in pp collisions as [12], 



dN 
dy 



(b) 



n(b) 



dN n 



1 + S(A 1 / 3 + B 1 / 3 ) dy 



(8) 



pp/dy = 0.8 In ^fs is the hadron multiplicity for 



where dN } 

nucleon-nucleon collisions at mid-rapidity [14] . The value 
of dN/dy is 725 at b ~ 3.2 fm, y 7 ! = 17.3 GeV and 
S = 0.09 (see Fig.2). The term 1 + 5(A^ 3 + B 1 ' 3 ) takes 
into account the energy degradation of the participat- 
ing nucleons in the nuclear environment. For isentropic 
expansion the initial temperature (Tj) is related to the 
hadron multiplicity as, 



Tfib) = 



2ir 4 



45£(3) 7ri%Ti4a fe dy 



d 4m 



(9) 



a fe = 7'" 2 ff/c/90 is determined by the statistical degener- 
acy (gk) of the system formed after the collision. In 
case of a deconfined initial state, gk = 37 considering a 
non-interacting gas of gluons and quarks of two flavours. 
For a hadron gas composed of w, p, u, i], a\ and nu- 
cleons the effective degeneracy has a value ~ 30 near 
the critical temperature where the hadronic masses go 
to zero according to the universal scaling scenario [5]. 
The increase in the degeneracy originates from the heav- 
ier hadrons going to a massless situation ( see also [15]). 
A value of g e ff ~ 30 can also be realized from the lat- 
tice data [16,17] near the phase transition point by us- 
ing the relation s/T 3 = 4ir 2 g e f f /90, where s is the en- 
tropy density. The latter is calculated using the relation, 
s = (e + P)/T and contains the effect of the in-medium 
masses of the constituent hadrons (for details see [9]). 



For Tj — 1 fm/c [18] and gk ~ 30, we get an initial tem- 
perature T t (b ~ 3.2 fm) ~ 200 MeV (see Fig. 2). Note 
that the variation of Ti(b) as a function of the impact 
parameter is very slow because T,(6) ~ [dN/dyty)] 1 / 3 . 
In an attempt to understand the lattice data in terms of 
the effective fields [19] the effective degrees of freedom 
in hadrons was limited to <~ 24 which coincides with the 
number of fundamental degrees of freedom for two quark 
flavors. 



4. SPACE-TIME EVOLUTION AND EQUATION 
OF STATE 

Next we study the sensitivity of the results on the 
space-time evolution. In order to do that we solve the 
(3+1) dimensional hydrodynamic equations with initial 
energy density [20], 



,(r-R A )/S + i 



(10) 



and initial velocity profile, v r = vq (r/RA) a which has 
been successfully used to study transverse momentum 
spectra of hadrons [21,22] and photons [23-25]. For our 
numerical calculations we choose a = 1 and Vq = at 
the initial time. 

The equation of state (EOS) used here to solve the 
hydrodynamic equations is evaluated as in Ref. [9]. The 
temperature dependence of the mass enters into the EOS 
through the effective statistical degeneracy {geff)- In 
Fig. 3 the temperature variation of the effective degen- 
eracy obtained from different models for the hadronic 
interactions (see [9] for details) is compared with the 
lattice QCD calculations [26]. The co-efficient of e/T 4 
differs from the effective degeneracy by a factor of 7r 2 /30. 
The lattice data seems to be well reproduced by the 
universal scaling scenario of mass reduction given by 
eq. 1. In this case the velocity of sound (c s ) is given by 
c- 2 = Tds/sdT = [(T/g eff )(dg eff /dT) + 3} [27]. Clearly, 
the value of c s is less than its value corresponding to 
ideal fluid case ( c l s deaL = This affects the space 

time evolution of the system non-trivially. In solving the 
hydrodynamic equations we have used the equation of 
state which contains the in-medium shift of the hadronic 
spectral functions as described above. 

The integration over the space time history has also 
been performed by taking the temperature profile from 
the transport model where the temperature varies with 
time as follows: 



T(t) = (T t -T 00 )e- t ' T + 



(11) 



The cooling law given above is the parametrization of 
the results of [28], used in several articles in the litera- 
ture [29] (see also the review [8] and references therein). 
The calculation is performed with the following values of 
parameters: T t = 200 MeV, = 120 MeV, r = 8 fm/c. 
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FIG. 3. The energy density e in the unit of T for the 
equation of state for various models of hadronic interactions 
is plotted as function of temperature (T) in the unit of the 
critical temperature, T c . The filled circle denotes the lattice 
results [26]. The universal scaling scenario seems to describe 
the lattice results quite well near the T c . For T > T c the bag 
model equation of state has been used. 



5. HADRONIC SPECTRA AND FREEZE-OUT 
TEMPERATURE 

Having fixed the initial temperature and the equation 
of state we now need to fix another important parameter, 
the freeze-out temperature Tp where the hydrodynamic 
evolution should terminate. At this stage the mean free 
path of the constituents begin to exceed the size of the 
system and particles start free streaming to the detector . 
For this purpose we will consider the mr(= y/p^+rn?) 
spectra [30,31] of pions and protons measured in Pb + Pb 
collisions at SPS energies by the NA49 collaboration [4]. 
We will assume that pions and nucleons are in thermal 
equilibrium throughout the evolution until they freeze- 
out at a common temperature Tp. Negative hadrons 
and positive minus the negative are treated as pions and 
protons respectively. The 7r~ and K~ from the decays 
p and (f> are about 5% of the direct pions at px ~ 500 
MeV and are therefore neglected here. In fig. 4 NA49 
pion spectra is compared with the hadronic initial states 
with Ti ~ 200 MeV and Tp = 100, 120 and 140 MeV. We 
find a reasonable agreement for Tp — 120 MeV as far as 
the slope is concerned. Calculations are done with all the 
scenarios I, II and III. Due to reasons explained in section 
2, the results with scenario II are indistinguishable from 
that of scenario III. We arrive at similar conclusions from 
the analysis of the proton spectra as shown in fig. 5. 
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FIG. 4. The mr distribution of pions for Pb + Pb collisions 
at CERN SPS energies. 
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FIG. 5. The mr distribution of protons for Pb + Pb colli- 
sions at CERN SPS energies. 



6. PHOTON SPECTRA 

We start with the single photon data of WA98 collab- 
oration which has initiated considerable theoretical ac- 
tivities [23-25,32-34]. For the analysis of WA98 photon 
spectra, we consider prompt photons resulting from the 
hard collisions of the partons in the colliding nuclei and 
the thermal photons from a hot hadronic gas with pos- 
sible in-medium modifications. Photons from hadronic 
decays (ir° — > 77, 77—8-77 etc.) are already subtracted 
from the data [2] and hence we need not consider them 
here. We begin our discussions with the prompt photon. 
The prompt photon yield for nucleus- nucleus collision is 
given by, 



d d p (7i n d 6 p 



(12) 



where n(b) is the average number of nucleon-nucleon col- 
lisions at an impact parameter 6, and Gi n (=30 mb, 
taken from Ref. [2]) is the p — p inelastic cross sec- 
tion. The prompt photon contributions have been evalu- 
ated with possible intrinsic transverse motion of the par- 
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tons [33,35] inside the nuclcon and multiplied by a K- 
factor <~ 2, to account for the higher order effects. The 
CTEQ(5M) parton distributions [36] have been used for 
evaluating hard photons. The relevant value of y/s, en- 
ergy in the centre of mass for WA98 experiment is 17.3 
GeV. No experimental data on hard photons exist at 
this energy. Therefore, the "data" at y/s — 17.3 GeV 
is obtained from the data at yfs — 19.4 GeV of the 
E704 collaboration [37] by using the scaling relation: 
Ed<r/d 3 p 7 \ hl+h2 ^ c+1 = f(x T = 2p T / v /~s)/s 2 , for the 
hadronic process, hi + h 2 — > C + 7 [38]. This scaling is 
valid in the naive parton model. However, such scaling 
may be spoiled in perturbative QCD due to the reasons, 
among others, the momentum dependence of the strong 
coupling, a s and from the scaling violation of structure 
functions, resulting in faster decrease of the cross section 
than 1/s 2 [35]. Therefore, the data at V« = 17-3 GeV 
obtained by using the above scaling gives a conservative 
estimate of the prompt photon contributions. We have 
seen that the effects of the nuclear shadowing in par- 
ton distributions on the prompt photons are negligibly 
small at SPS but it is important at RHIC and LHC en- 
ergies [39]. This is because the value of ie(<~ 2px/^/s) at 
SPS is not small enough for the shadowing effects to be 
important. 

To evaluate the photon emission rate from a hadronic 
gas we model the system as consisting of n, p, u r\ and 
a\ [40,41]. The relevant vertices for the reactions irir — > 
p 7 and 7r p — > 7r 7 and the decay p — > tt tt 7 are obtained 
from the following Lagrangian: 

C = -g P ™W ■ (*? x M 7f) - eJM M + (p„ x p v ) 3 , 

(13) 

where F^ v = d p A v — d v A pi is the field tensor for electro- 
magnetic field and J p is the hadronic part of the electro- 
magnetic current given by 

J« = {Pu x ^) 3 + (t? x (c^tt + g pnn n x ^)) 3 , (14) 

with q^ v = d p p v - d v p p - g p ^{p p x p v ). tt, p^ and A* 
represent the tt, p and photon fields respectively and the 
arrows represent vectors in isospin space. g pnn denotes 
the coupling strength of the pirir vertex, fixed from the 
observed decay width p — > nn. We have also considered 
the photon production due to the reactions tt rj — > 777, 
tt tt — > 77 7 and the decay to — > tt 7 using the interaction 
given in [42] and vector meson dominance [43]: 

+ £ ^V^A, (15) 
v= P ^ 9v 

where e pua fj is the totally anti-symmetric Levi-Civita 
tensor. The invariant amplitudes for all the reactions 
are given in Ref. [41]. The values of g ppri and g wp7r are 



fixed from the observed decays, p — > 777 and oj — > 717 
respectively. The constant gy is determined from the de- 
cays, V — > e + e~ . Photon production due to the process 
7r p — ► a\ — » 7T7 is also taken into consideration (for 
details of interaction vertices see [9]). 

In Fig. 6 the static (fixed temperature) photon spec- 
tra is shown for T = 180 MeV. The solid line shows 
the enhanced yield with in-medium masses and vac- 
uum widths compared to the yield obtained with vac- 
uum masses and widths (dotted line). The long- 
dashed line showing the result with in-medium width 
and vacuum masses does not differ substantially from 
the dotted line, because we find that the effects of the 
change in the decay width both on the phase space 
and the production cross section are small. We ob- 
serve that the contributions from the decays of bary- 
onic resonances (V(1520), V(1535), V(1440), A(1232), 
and A(1620)) are also small (filled circle). The values 
of the decay widths, R — > N 7, where R and N denote 
the baryonic resonances and the nucleon respectively, are 
taken from the particle data book [44]. It is observed 
that the photon yield with universal mass variation sce- 
narios, I is almost an order magnitude larger than the 
case II with large broadening of the p. Reduction in the 
hadronic masses causes an enhancement in the thermal 
distribution of mesons because of which the rate of ther- 
mal emission of photons increases. The scenario with 
vacuum values of both masses and width (dotted line) 
does not show any appreciable difference from (II). 
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FIG. 6. The photon production rate as function of its en- 
ergy at a temperature, T = 180 MeV. Solid (long dashed) line 
indicates rate for scenario I (scenario II). Dotted line repre- 
sents the spectra without any medium effects. Filled circles 
indicate photons originating from baryonic resonance decays. 
Short dashed (dot-dashed) line indicates the photon emission 
rate from a thermalized two-flavor quark gluon plasma for 
a s = 0.08(0.2). 
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For the sake of comparison we have also shown the 
photon emission rate from a thermalized QGP. The pho- 
ton spectra from QGP includes Compton, annihilation, 
bremsstrahlung and annihilation with scattering pro- 
cesses [40,45-47] (see also [48]). Here we have mainly 
focused on photon production from hot hadronic mat- 
ter and hence do not consider emission rate from QGP 
beyond two loops [49]. From the lattice results [50] the 
value of a s ~ 0.3(y 2 = 4) at T = 180 MeV. At T = 180 
MeV the emission rates from the QGP (dot-dashed) and 
hot hadronic matter (scenario II and III) are comparable 
for pt > 1.5 GeV. The photon yield from hadronic matter 
for scenario I is substantially larger than that from QGP. 
However, it should be mentioned here that the photon 
emission rate from QGP is evaluated in Refs. [40,45,46] 
by using hard thermal loop (HTL) approximations, which 
is valid for the value of the color charge, g << 1. The 
validity of the HTL approximation to higher values of g 
is doubtful. In fact, such an extrapolation in the value 
of g (from g « 1 to g ^2) introduces an uncertainty 
by a factor of <~ 5. Due to these uncertainties we must 
be careful in making any firm conclusion from the re- 
sults obtained on the basis of the HTL approximations. 
For a s = 0.08(g = 1) the emission rate from QGP (short- 
dashed line) is much smaller than the rates from hadronic 
matter. 

In Fig. 7 the pt distribution of 
photons (prompt+thcrmal) is compared with the WA98 
data. Within the framework of the transport model the 
data is well reproduced when the hadronic masses are 
allowed to vary according to the eqs. 1 (long-dash line, 
scenario I). However, when scenario (II) is considered for 
thermal photons (dash dotted line) the experimentally 
observed "excess" photon in the region 1.5 < pr (GeV) 
< 2.5 is not reproduced. The dotted line indicates results 
with vacuum masses and widths. 
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FIG. 7. Total (prompt+thermal) photon yield in Pb + Pb 
collisions at 158 A GeV at CERN-SPS. The theoretical cal- 
culations contain hard QCD and thermal photons. The sys- 
tem is formed in the hadronic phase with initial temperature 
Ti = 200 MeV; 'trans' indicates the results for the cooling 
law 11. 

The photon yield for the scenario (I) shown in Fig. 7 
(solid line) is obtained by folding the static rate with the 
space time history of the system governed by relativistic 
hydrodynamics from the initial to the freeze-out state 
(freeze-out temperature ~ 120 MeV); the data is well 
reproduced in this case also. 

The agreement between the px spectra of photon ob- 
tained with two different types of space time evolution 
scenarios (transport model and hydrodynamics) can be 
explained as follows. We find that the variation of tem- 
perature with time (cooling law) in cq. 1 1 is slower than 
the one obtained by solving hydrodynamic equations. As 
a consequence the thermal system has a longer life time 
than the former case, allowing the system to emit pho- 
tons for a longer time. In case of hydrodynamics this is 
compensated by the transverse kick due to radial velocity 
of the expanding matter which shifts some of the photon 
multiplicity from lower values of transverse momentum 
towards the higher region. 



7. DILEPTON SPECTRA 

Now we study the invariant mass distribution of lep- 
ton pairs measured by CERES/NA45 collaboration [3] 
in Pb + Au collisions. We consider the thermal dilepton 
production in the hadronic medium due to the process 
7r + 7r~ — » e + e~, known to be the most dominant source 
of dilepton production from hadronic matter. In the ther- 
mal system the width of the to meson can be large due to 
various reactions occurring in the thermal bath, the most 
dominant process, among others, is lutt <-> 7T7T [9,51]. As 
a consequence of this, the life time of the u> meson could 
become smaller compared to the life time of the hadronic 
system, enabling it to decay in the interior of the sys- 
tem. In view of this the in-medium decay uj — > e + e~ 
is also taken into account. The required interaction ver- 
tices have been obtained from eqs. (13), (14) and (15). In 
Fig. 8 the experimental data is compared with the the- 
oretical results for dN c h/drj = 270. The effective masses 
and widths of p and uj mesons appearing in the processes 
-kit — > p — > e+e~ and uj — > e + e~ respectively are taken 
from Eq. 1 (for details on the the thermal emission rate 
of lepton pairs see Ref. [9]). Contributions from hadronic 
decays at freeze-out (background) is taken from the sec- 
ond of Ref. [3] (dotted line). The observed enhancement 
of the dilepton yield around M ~ 0.3 — 0.6 GeV can be 
reproduced by the scenarios (I) (long-dashed) and (II) 
(dot-dashed) when the cooling law is taken as eq. 11. 
The data is also well reproduced when hydrodynamics is 
used to describe the space time evolution of the system 
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for the scenario I (solid line). However, with vacuum 
properties of the vector mesons the low mass enhance- 
ment can not be reproduced (short-dashed), indicating 
the change of vector meson properties in the medium. It 
is interesting to recall here that in a recent experiment 
Ozawa et al [52] has observed a significant enhancement 
in the dilepton yield in p + Cu collisions as compared to 
p + C collisions below the ui peak. 




M (GeV) 

FIG. 8. Dilepton spectra for (N c h}=270 for different sce- 
narios as indicated in the text. 



8. SUMMARY AND DISCUSSIONS 

Let us now summarise the results. We reiterate that 
our main contention in this work is to comment on the 
nature of medium effects vis-a-vis mass shift as opposed 
to broadening of the spectral density of hadrons. Such a 
study is only possible through the electromagnetic probes 
since they are sensitive to the evolution of the spectral 
function of vector mesons. The freeze-out conditions 
within the model used for evaluating the photon and 
dilepton spectra have been estimated through the NA49 
pion and proton spectra. Firstly, we find that the WA98 
photon transverse momentum spectra and the CERES 
dilepton low invariant mass spectra is well explained if 
we assume a chirally restored initial phase at a tempera- 
ture <~ 200 MeV where the masses of the hadrons tend to 
zero and grow to their vacuum values with the decrease 
of temperature. In fact, close to the critical temperature 
for chiral symmetry restoration (in our case ~ Tj = 200 
MeV) the description of the system either in terms of 
partonic or hadronic degrees of freedom become dual, al- 
though the non-perturbative effects may still be impor- 
tant. For example, the emission rate of dilcptons from 
pion annihilation will resemble that from qq annihilation, 
because of the complete extinction of the intermediary p 
at the critical temperature [8,9]. This indicates that the 
qq interaction in the vector channel has become weak, 
signalling the onset of chiral transition. Whether such 
a state is synonymous with deconfined matter (QGP) is 
still an unsettled issue. A similar value of the initial 
temperature is obtained in Refs. [24,25,27,53]. Based on 



the present and our earlier results [24], it is fair to say 
that a simple hadronic model with vacuum properties of 
hadrons is inadequate to explain the above experimen- 
tal data. Either a substantial change in the in-medium 
hadronic spectral function or the formation of the QGP 
is required. More importantly, we find that the pr dis- 
tribution of photons changes significantly with a reduced 
mass scenario and is almost unaffected by the broaden- 
ing of the vector meson spectral function in the medium. 
The invariant mass distribution of the lepton pairs on 
the other hand can be explained with both a reduction 
in the mass as well as with an enhanced width of the 
vector mesons. Thus by looking only at the dilepton 
spectra it is difficult to differentiate the above scenarios; 
we need to analyse both the photon and dilepton spectra 
simultaneously. With better statistics we might succeed 
in ruling out or at least restricting one of the scenarios. 
We expect that such a situation should be realized at the 
Rclativistic Heavy Ion Collider (RHIC). 
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